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We present our recent studies on hyperon-nucleon (YN) interactions in the strangeness 5 = — 1 that 
pA,lPp and H+n, by extracting corresponding potentials through Nambu-Bethe-Salpeter wave 
functions. We calculate AA^ and EA^ potentials in the isospin / = 3/2 channel, using the Nf — 
2+1 gauge configurations generated by PACS-CS collaboration and employing an improved 
method to obtain potentials in lattice QCD simulations. For the '^o channel, the central EA^(/ = 
3/2,^ So) potential and the central AN{^Sq) potential are found to be very similar. In the spin 
triplet i^Si —^Di) channels, the central AN{^Si —^Di) potential is attractive while the central 
EA^(/ = 3/2,^ ^i —^Di) potentials is repulsive. Tensor potentials, on the other hand, are rather 
weak in the diagonal part of both AN and EA^(/ = 3/2) systems. 
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1. Introduction 

The A-nucleon {AN) and the Z-nucleon {l^N) interactions are one of the basic inputs to study 
the hypemuclear systems, in which hyperons (or strange quarks) are embedded in normal nu- 
clear systems as "impurities" [^. For example, from studies of few -body systems for 5-shell A- 
hypernuclei, it has been pointed out that a coupled-channel AN — LN interaction plays a significant 
role to make a hypemucleus being bounded. [^]. On the other hand, no Z-hypernuclei have been 
observed except for four-body Z-hypemucleus (^He), and furthermore a recent experimental study 
suggests that Z-nucleus interaction is repulsive. Such informations are useful to study properties 
of hyperonic matters inside the neutron stars [||], though a hyperonic equations of state (EOS) em- 
ployed in such a study may contradict a recent observation of a massive neutron star heavier than 
2Mq I^. Despite their importance, present phenomenological AN and JIN interactions have still 
large uncertainties since direct AN and LN scattering experiments are either difficult or impossible 
due to the short life-time of hyperons. 

Recently, a new approach has been proposed lattice QCD to study not only the NN interaction 
but also various baryonic interactions including the H-dibaryon system [|5], ^, ^, |8[ ^ 10, |n|], and 



the method is recently extended to systems in inelastic channels [ p^ |T3|]. This method has been 



applied also to study various interactions other than baryon-baryon interactions [ ]14| , |15| ]. In this 
approach, the interparticle potential can be directly extracted in lattice QCD through the Nambu- 
Bethe-Salpeter (NBS) wave function, and observables such as the phase shift and the binding en- 
ergy can be calculated by using the resultant potentials. The purpose of this report is to present 
our recent calculations of the AN potentials as well as the LN{I = 3/2) potentials using full QCD 
gauge configurations. Several earlier results had already been reported at LATTICE 2008 [|T^] and 
LATTICE 2009 [jT/l]. This report is the latest version of those reports, which includes several new 



results: (i) l.N{I = 3/2) potentials are studied, (ii) An improved method is employed to extract 



potentials more precisely [18] 



2. Improved extraction of potentials 

In the HAL QCD scheme, the non-local but energy-independent potential is defined by the 
Schrodinger equation as 

(V2 + ^2) <^{r)=2iijd\'U{r,r')<^{¥), (2.1) 

where 0(r) is the equal-time NBS wave function of two-baryon system (61,62). M ="JBi'WB2/("^Bi + 
m^j) and k^ are the reduced mass of the (81,82) system and the square of asymptotic momentum 
in the center-of-mass frame, respectively. In practice, the nonlocal potential is expanded in terms 
of the velocity(derivative) as [p^], l]{r,r') = Vg[B,(r, V)5(r — r'). The potential V may have an 
antisymmetric spin-orbit term when two baryons are not identical. For example, the AA^ potential 
is given by 

Vkn = VQ{r)+Vc{r){aK-aN)+VT{r)Sn + VLs{r){L-S+)+VALs{r){L-S-) + 0{W^), (2.2) 

where 512 = 3(aA-n)(aA?-H) — df^-OM is the tensor operator with n = r/\r\,S± = {GNrtdfO/l are 
symmetric (+) and antisymmetric (— ) spin operators, L = —if x V is the orbital angular momentum 
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operator. In the velocity expansion, Vb.fj.r are the leading order (LO) potentials, denoted by Vlo. 
while VisALS are of the next-to-leading-order (NLO). In this report, we consider the LO potentials 
only. 

In lattice QCD simulations we first calculate the normalized four-point correlation function 
defined by 



<r(^-^-^o)=i 



{J.M) , 



Bi,«(X + r,0B2,/3(^,0XB4 (^o) 



0Wexp{-(mBi+mB2 )(?-?())}, 



(2.3) 
where the summation over X selects states with zero total momentum. The B\^a{x) and B2^j}{y) 
denote the interpolating fields of the baryons such as 



P = £abc (UaCYsdh) Uc, U = Sabc {UaCy^db) 4, (2-4) 

S^ = -£ahc{UaCj5Sb)Uc, iP = —7=£abc{{daCj5Sb)Uc + {UaCj5Sb)dc], (2.5) 



A = —7=eabc{{daCy5Sb) Uc + (SaCysUb) dc - 2 {UaCjsdb) Sc} ■ 



(2.6) 



rilM). 



,(J.M)- 



and ^g ^ (?o) = La'j3' ^a'R' ^3,0;' (?o)^4./5' (^o) is a source operator which creates B3B4 states with 
the total angular momentum J,M. This normalized four-point function can be expressed as 



JJ,M) 



£,A e^^^""'""! -"'fi2)('-'o), (2.7) 



where £"„ (!£■„)) is the eigen-energy (eigen-state) of the six-quark system with the particular quan- 
tum number (i.e., J^,M, strangeness S and isospin /), and A„ = La'S'^a'R' (^n|^4,i3'^3,a'|0)- 



Since E„ - mb, - vub^ = k^/{2}x) + 0{k'^), we have[|r8|] 

R{f,t) = ld\'U{r,r')R{r',t) + 0{k^)=VLoir)R{r,t) + ---, (2.8) 



IjX dt 



where t — to should be moderately large so that states with large k^ and states with more than 2 
particles are suppressed. 

For the spin singlet state, we extract the central potential as Vc(r;7 = 0) = (^ — ^)R/R. For 
the spin triplet state, the wave function is decomposed into the S- and the D-wave components as 



I Rapin 'Si) = ^Rap(r;J=l) = JrA^MeO^Rap{rJ=l), 
\ Rap{r, 'Di) = £2Rap{7;J = 1) = (1 - ^)Rap{r,J = 1). 

Therefore, the Schrodinger equation with the LO potentials for the spin triplet state becomes 



(2.9) 



^ ( ^<y-:p + VQ{r)+Vcy{r){di^-dN) + VT{r)Sn\R{7,t-to) 



3^ d 

(2.10) 



from which the central and the tensor potentials, Vc{r,J = 0) = Vo{r) — 3V(j{r) for 7 = 0, Vc{r;J 
1) = V'o(r) + Va{r), and Vrir) for / = 1, can be determined. 
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Figure 1: Left: The central potential (red) in the 5o channel of the AA^ system in 2 + 1 flavor QCD as a 
function of r, together with the ^ part (green), 4- part (blue). Right; The central potential in the ' ^o channel 
of the ILNil = 3/2) system as a function of r. Symbols are same as the left figure. 



3. Numerical simulations 



In our calculations for AN and l.N{I = 3/2) potentials, we employ 2-1-1 flavor full QCD gauge 



configurations generated by PACS-CS collaboration [23| with the RG-improved Iwasaki gauge 
action and the nonperturbatively O (a) -improved Wilson quark action at j3 = 1.9 on a L^ x T = 
32^ X 64 lattice. The lattice spacing at the physical quark masses is a = 0.0907(13) fm [^], 
thus the spatial volume is (2.90fm)^. We have chosen one set of hopping parameters {Kud, k:^) = 
(0.13700,0.13640) for light and the strange quarks, which correspond to {mj^niK) = (699,787) 
MeV. The wall source with the Coulomb gauge fixing is placed at the time-shce Iq, and the Dirichlet 
boundary condition is imposed in the temporal direction at the time-slice t — to = Nt/2. A total 
number of gauge configurations we have used is 399, and we put a source at to = Sn with n = 
0, 1,2, • • • ,8 on each configuration to increase the statistics. 



4. Results 

4.1 AA^ and LN{I = 3/2) potentials in ^Sq channel 

The AN (left panel) and the rA'^(/ = 3/2) (right panel) potentials in the '^o channel are shown 
in Figure ^ where the laplacian part is represented by green, the time-derivative part by blue and the 
total potential by red, respectively. Note that the time derivative term gives a significant attraction 
at medium and longer distances. 

In the 2-1-1 flavor QCD, while the LN (I = 3/2) potential stiU belongs directly to the 27(7 = 
3/2) representation thanks to the isospin(/) symmetry, an energy eigenstate of a AN system in 
the ^^0 channel is a mixture of 27(7 = 1/2) and 8.v in the flavor representation, so that these two 
potentials are not necessarily mutually in agreement. As seen from Figure |I], however, these two 
potentials look very similar to each other, since the flavor symmetry breaking is still small in 
the present 2 + 1 flavor QCD calculation, where the baryon masses are given by {mN,mA,mz) = 
(1.574(3), 1.635(3), 1.650(3)) GeV. 
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Figure 2: Left: The central potential in the ^^i —^D\ channel of the AA^ system as a function of r. Symbols 
are same as the Fig. |l|. Right: The tensor potential in the ^S\ —^ Di channel of the AA^ system as a function 
of r. Symbols are same as Fig. nl 



4.2 AN interaction in the ^^i —^Dy channel 

Figure ^ shows the central potential (left panel) and the tensor potential (right panel) of the AA^ 
system in the ^^i — ^Di channel, whose eigenstate is a mixture of 10 and 8^. The time-derivative 
part (blue) of the central potential has attractive contributions, so that the attraction of the total 
potential (red) is enhanced at long distance. The attractive well at the distance r ss 0.6 fm is deeper 
than that of the AN central potential in the ^5o channel. The time derivative part of the tensor 
potential, on the other hand, gives a negligible contribution. Furthermore, the tensor potential itself 



(red) is weaker than the tensor potential of the NN system| 18]. 



4.3 ILN interaction in the ^Si —^D\ channel 

Figure ^ shows the central potential (left panel) and the tensor potential (right panel) of the 
'LN{I = 3/2) system in the ^Si —^Di channel. Due to the isospin symmetry, this channel belongs 
solely to the flavor 10 representation without mixing with 10 or 8^ 

As seen form an enlargement in the inset of the left panel, there is no clear attractive well in 
the central potential (red). An attractive well seen in the laplacian part (green) at medium distances 
is cancelled by a repulsive contribution in the time-derivative part (blue) at the whole range. This 
repulsive nature of the rA'^(/ = 3/2,^ Si —^Di) central potential is consistent with the prediction 
from the naive quark model. |^^ The tensor force is a little stronger that that of the AN system but 
is still weaker in magnitude than that of the NN system. The time-derivative part gives a negligible 
contribution. 



5. Summary 

We have calculated the central and the tensor parts of the AA'^ and rA^(/ = 3/2) potentials using 
the improved method. Two potentials in the ^Sq channel, one is AA'^ and another is £A'^(/ = 3/2), 
are very similar to each other, since both sectors have a common 27 representation in the flavor 
SU(3) and the SU(3) breaking effect is still small in the present calculation. Both potentials have 
an attractive well, so that they give attractive interactions at low energy. 
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Figure 3: Left: The central potential in the ^^i —^ D\ channel of the EA^ system as a function of r. Right: 
The tensor potential in the ^S\ —^Di channel of the 'LN system as a function of r. 



In the ^^i —^Dy channel, while the central AN potential shows attraction at low energy with 
an attractive well, the central 'LN{I = 3/2) potential is repulsive at all distances. Tensor potentials 
for both AA^ and rA'^(/ = 3/2) systems in this work are rather weak. 

In this report, we have extracted single-channel potentials only without coupled-channel analy- 
sis between AN and 1.N . though the coupled-channel potentials can be also extracted in principle! 12, 
13]. The author had postponed such a coupled-channel analysis to future publications, since he had 
to spend his time to recover in various parts of his life and research environment from the massive 
earthquake in East of Japan on March 11th, 2011. Further analyses including coupled-channel 
potentials are now in progress. 
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